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ABSTRACT 

We have performed Smoothed Particle Hydrodynamic (SPH) simulations to study 
the response of molecular clouds in the Galactic disk to a rotating bar and their 
subsequent evolution in the Galactic Center (GC) region. The Galactic potential in our 
models is contributed by three axisymmetric components (massive halo, exponential 
disk, compact bulge) and a non-axisymmetric bar. These components are assumed to 
be invariant in time in the frame corotating with the bar. Some noticeable features 
such as an elliptical outer ring, spiral arms, a gas-depletion region, and a central 
concentration have been developed due to the influence of the bar. The rotating bar 
induces non-circular motions of the SPH particles, but hydrodynamic collisions tend 
to suppress the random components of the velocity. The velocity field of the SPH 
particles is consistent with the kinematics of molecular clouds observed in HCN (1 — 0) 
transition; these clouds are thought to be very dense clouds. However, the I — v diagram 
of the clouds traced by CO is quite different from that of our SPH simulation, being 
more similar to that obtained from simulations using collisionless particles. The I — v 
diagram of a mixture of collisional and collisionless particles gives better reproduction 
of the kinematic structures of the GC clouds observed in the CO line. The fact that 
the kinematics of HCN clouds can be reproduced by the SPH particles suggests that 
the dense clouds in the GC are formed via cloud collisions induced by rotating bar. 

Subject headings: The Galaxy: center - clouds: structures : SPH: gas dynamics 
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1. INTRODUCTION 

The properties of molecular clouds in the Galactic Center (GC) are known to be influenced by 
special environmental effects such as strong tidal fields, strong mass concentration and energetic 
activities. There have been many observations to study the structure and kinematics of such 
GC molecular clouds. Since different molecular lines trace clouds of different densities, the 
comparison of several different observations should be very useful in understanding the nature 
of the clouds in the GC. For example, CO (Bitran 1987) and HCN (Lee 1996) data show quite 
distinct features. The HCN data give informations on the distributions and dynamical structures 
of "dense" 10 5 cm~ 3 ) molecular clouds, whereas the CO line traces "less dense" (<^ 10 3 cm -3 ) 
clouds. In CO (1 — 0) observations, it is difficult to see structures with velocities |Vl5_r| ~ 10 
km s _1 toward the GC region because of the heavy contamination by the disk components which 
have small line-of-sight velocities. But in the HCN (1 — 0) line, there is no contamination in this 
velocity range because there is very little HCN emission from the local disk. On the other hand, 
with HCN data alone, we lose important information about certain interesting features of the 
low-density clouds in the GC such as the high-velocity (|Vls.r| ~ 170 — 260 km s^ 1 ) components. 
Consequently, the I — v maps in CO and HCN are quite different. For example, the CO data show 
a clear parallelogram structure (Binney et al. 1991) which is not seen in the HCN maps. Thus, 
surveys with high-density and low-density tracers serve us with complementary information on the 
distribution and kinematics of GC molecular clouds. Other molecular lines such as CS also trace 
the dense clouds. However, we consider only HCN observations for dense components because the 
CS surveys by Bally et al. (1987) and Tsuboi et al. (1989) cover a smaller area than the HCN 
survey. We simply note that the CS and HCN data are very similar in overlapping areas. 

One of the interesting features of the CO and HCN maps includes the existence of clouds 
with "forbidden velocities" (velocities that differ significantly from what would be observed along 
a given line-of-sight if the clouds are in purely circular orbits). It would be difficult to explain 
these non-circular velocity components if the Galactic potential is axisymmetric, unless there is a 
mechanism for generating non-circular motion such as explosive phenomena. 
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There are direct evidences for a bar in the inner Galaxy (e.g, non-axisymmetric image of the 
GC region seen from COBE satellite; Blitz et al. 1993) as well as somewhat indirect evidences 
(e.g, parallelogram structure of CO in the I — v map) deduced from the kinematic behavior of 
molecular clouds [see Gerhard (1995) for a review]. Thus it is desirable to examine the effects of 
the Galactic bar on the distribution and kinematics of molecular clouds. 

There have been many hydrodynamic simulations for the study of interstellar medium in the 
Galaxy: sticky particle simulations by Schwarz (1981, 1984) and Habe & Ikeuchi (1985), smooth 
fluids simulations by Sanders & Huntley (1976), van Albada (1985), Mulder & Liem (1986), and 
Athanassoula (1992), and Smoothed Particle Hydrodynamic (SPH) simulations by Wada & Habe 
(1992), Friedly & Benz (1992), and Fux & Friedly (1996). These simulations have been successful 
in reproducing gross structures of the interstellar medium and gas motions on global (a few kpc) 
scales, but there has not been much effort in the detailed study of the structures within ~1 kpc 
from the GC. In the present study, we concentrate on the kinematic structures in such small 
region of the GC. We pay special attention to the different structures revealed by CO and HCN 
observations. 

Recently, Binney et al. (1991) have shown that single-particle orbits in a steady non- 
axisymmetric potential should provide a reliable guide for the interpretation of kinematic behavior 
of molecular clouds in the GC region. They fitted successfully the large non-circular forbidden 
components with the CO parallelogram in the I — v map using xl orbits. They also showed that 
locations of the giant molecular clouds at the GC, such as the Sgr B and Sgr C, could be explained 
by the x2 orbits. However, more realistic calculations including the non-axisymmetric potential 
as well as the hydrodynamic effects are highly desirable to confirm some of the conclusions 
drawn from orbit calculations. For such purposes, we carry out numerical simulations using the 
SPH method for the response of molecular clouds in the Galaxy with a bar. We then make 
detailed comparisons between the observation and simulations to understand the distribution and 
kinematics of the GC molecular clouds. 

In § 2 we shortly introduce the observational data which will be compared with our simulation 
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results. In § 3 we describe the model Galaxy, model parameters, and the gaseous disk. The 
results from the numerical calculations are compared with ones observed in the CO and HCN 
molecular lines in § 4. In § 5 a run of mixture of collisional and collisionless particles, the gas infall 
phenomenon driven by the perturbation of bar potential, and its related structures are discussed. 
Finally, we summarize the results in § 6 

2. OBSERVATIONAL DATA 

The kinematic behavior of the clouds can be studied using I — v maps of molecular emission 
lines. In Fig. |l|, we have shown CO and HCN I — v maps. The CO I — v map (Fig. [l]a; Bitran's 
data kindly provided by T.Dame) shows a clean parallelogram structure with a maximum velocity 
of around 300 km s _1 . The emission resulting from clouds in the disk component (from solar 
neighborhood to near the center) is seen as a zero velocity 'strip' in this map. There are relatively 
empty regions near I ~ ±(2 — 4°). The emission is stronger at positive I than negative I. The HCN 
map (Fig. ||b) is based on the large scale survey carried out by Lee (1996) at the Taeduk Radio 
Astronomy Observatory. The HCN emission is generally much weaker than the CO and it is very 
strongly skewed toward the positive I region. The HCN map does not show the parallelogram 
structure seen in the CO map. The forbidden velocity components are much less prominent in the 
HCN map. In the following numerical study, we try to understand these differences between the 
CO and HCN data. 

3. MODELS 

There is no single numerical method that can describe the evolution of the galaxy since there 
are many components that react differently. We have adopted the SPH method to study the 
dynamics of molecular clouds in the Galaxy because of its versatility. The gravitational potential 
of the Galaxy is generated by the stars, gas, and dark matter. The Galactic potential is assumed 
to be independent of time (in the frame corotating with the bar) for simplicity. The molecular 
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clouds represented by SPH particles change with time, and so does the potential generated by the 
clouds. However, the total mass contained in this component is a small fraction of the Galactic 
mass and thus can be ignored in the first approximation. Simulations with self-gravity are highly 
desirable, but they require much better numerical accuracy and computational resources than the 
simple calculations reported in this paper. 



3.1. The Model Galaxy 

The model Galaxy is assumed to be composed of four components: halo, disk, bulge, and bar. 
We have ignored the central black hole (Eckart & Genzel 1996) because its influence is limited 
within about a parsec from the center. Our computational resolution is of order of 100 pc and the 
central black hole would have no significant effect in our simulations. 

Let's denote r as the distance from the GC to any location in the Galaxy and R as the 
distance on the disk. The halo component, which gives rise to the flat rotation curve at large radii, 
is assumed to have a logarithmic potential, 

Hr) halo = -v\ ln(R 2 c + r 2 ) + const, (1) 

where R c is the halo core radius, and vq is the constant rotation velocity at large r. For the disk 
component, we adopt Freeman's (1970) exponential disk which has the following form: 



HR)disk = vtGSq^ 



2R d ) k ° GtfJ h Gi?J Ki G;?J 



(2) 



Here Rd is the disk scale length, So is the central surface density, and Jo, Kq, I\, and K\ are 
modified Bessel functions. The total disk mass is simply Mdisk = 27rSoi?^. As for the bulge 
component, we have assumed the Plummer model to allow a steep outward velocity increase from 
the GC (Binney & Tremaine 1987); 

where r c is a parameter which controls the size of bulge. The bulges of real galaxies can be better 
represented by de Vaucouleurs' i? 1 / 4 -law, but we have chosen the Plummer model for simplicity. 
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Again, our results are not very sensitive to the specific choice of &buige as long as the typical radius 
where bulge dominates the potential is small. Our models have r c ~ 230 — 240 pc. 

The bar is a triaxial component in three dimensions. However, our simulation is restricted to 
the two-dimensional disk. Thus we have adopted the following form of the gravitational potential 
proposed by Long & Murali (1992), 

GM bar x-a + T_ 

HR)bar = ^ log — — ), 

2a x + a + 1 + 

where T± = y/[(a =L x) 2 + y 2 + b 2 ], a and b are the major and minor axes, respectively. In 
the above notation, x and y are the coordinates parallel to the long axis and the short axis, 
respectively. We have considered four models with bars of different axial ratios of a:b=2:l, 2.5:1, 
3:1, and 4:1. 



3.2. Model Parameters 

In our numerical calculations all physical quantities are given with dimensionless units. The 
unit of length in our calculation is chosen to be R sc =10 kpc. The total mass of the Galaxy within 
10 kpc is assumed to be M w = 1.24 x 10 11 M (Caldwell & Ostriker 1981). We have used the 
following units for time, velocity and angular velocity; 

time: t sc = \J ~ 424 x lo7 y ears > 

velocity: V sc = y w 230.7 km s _1 , 

angular velocity: U sc = w 23.07 km s _1 kpc -1 . 

The free parameters in defining the gravitational potential are chosen to satisfy the constraint 
that the rotation curve of the model Galaxy should match the observed one (Fig. |2|) . The observed 
rotation curve is taken from Clemens (1985). The mass of each component is assumed to be 
Mdisk ~ O.SljM&u/ge £3 0.07, Mb ar w 0.05, and the disk scale radius Rd is assumed to be 3.4 kpc. 
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We have considered four different values for b/a: 1/2, 1/2.5, 1/3 and 1/4. We have varied the axis 
ratios by changing the minor axis. The mass of each component has to be varied by small amount 
with axial ratio in order to best fit the observed rotation curve. The parameters for the models 
investigated in this study are summarized in Table 1. The final row in this table designated by 
'UB' denotes parameters of a model without a bar. 

The major axis length a of the bar, the bar pattern speed Qj, and the rotation speed v Q at the 
galactocentric distance of 10 kpc are fixed at 0.22 (Weiland et al. 1994; Dwek et al. 1995), 2.731 
(Gerhard 1995), and 0.92 (Caldwell & Ostriker 1981), respectively. The adoption of Q b = 2.731 
places the CR around at R =2.7 kpc, the ILR near R =1.2 kpc and the OLR close to R = 6 kpc 
(Fig- 1). 

3.3. The Gaseous Disk Model 

The extent of the exponential disk is much larger than the galactocentric distance of the Sun. 
However, we have distributed the SPH particles within 4 kpc in order to obtain better resolution 
with a finite number of particles. The number of particles for our simulations is around 8000. The 
initial resolution length was about 160 pc. 

We have used the SPH technique which is a Lagrangian method that uses a number of 
particles in order to solve the hydrodynamic equations. This is a versatile technique that can be 
applied to many three-dimensional problems. We simply refer readers to the review by Monaghan 
(1994) for detailed description of the SPH method. The code we used was adopted from J. 
Monaghan. 

In order to treat the hydrodynamics accurately, one needs to solve the the energy equation 
simultaneously. However, it is a very difficult task to include energy equation. The simplifying 
assumption to avoid such a complexity is to assume isothermality. The temperature T for the 
SPH particles can reach at high value when collision occur (which is highly supersonic in our 
cases). The resulting temperature is not a realistic gas temperature, but merely represents the 
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temperature just behind the shock. Even though the shock temperature can be fairly high, the 
cooling will be faster than any other dynamical processes. Therefore, we can assume that the 
shock is 'radiative' or 'isothermal'. 

In our numerical simulations, we have fixed the gas temperature at 10,000 K. Again, this 
temperature should be interpreted as the internal energy that is responsible to support the cloud 
against the gravitational collapse. Giant molecular clouds are known to support themselves mainly 
through turbulent motions. The velocity dispersion of about 10 km s~ 1 found in the molecular 
clouds roughly corresponds to temperature of 10,000 K. This assumption of constant temperature 
is based on the work by Hernquist (1989) who showed that the gas temperature can be maintained 
at ~ 10 K throughout the evolution due to an effective radiative cooling in the shocked gas. 
Clearly this choice is somewhat arbitrary, but many of the features of our numerical simulations 
are not sensitive to this temperature as long as the thermal energy is much smaller than the 
kinetic energy associated with the orbital motion. 

The initial randomly distributed velocities were given to the SPH particles so that the 
centrifugal accelerations of the particles were balanced with their gravitational accelerations. 

4. RESULTS 

4.1. The General Evolution of the Gas Distribution 

In order to make a test of our numerical code, we have first computed the UB model that 
does not have a bar, as a reference model. In this model there is little change in the particle 
distributions throughout the integrations because most of the particles keep circular orbits even 
at the end of simulations. 

In the models with the bar (B models), the orbits of the particles are much perturbed by the 
rotating bar potential and some noticeable features develop. In Fig. ^j, we show the distribution 
of SPH particles in reference frame at several different epochs for the B4-1 model. A two-armed 
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spiral pattern becomes visible at the early stages of the evolution (at t ~ 1) and remains weak 
to the later stages (t > 10). Other prominent features of the B models are the developments 
of non-circular orbits. The SPH particles near the CR experience large amount of velocity 
perturbations and eventually undergo dissipation of their orbital energy through hydrodynamic 
collisions. The inflow of mass due to this effect makes several interesting structures such as an 
elliptical outer ring at around 2 — 3 kpc, a prominent gas-depletion region between about 1 — 2 
kpc, and a central concentration of gas within about 1 kpc. 

4.2. Comparison with Observations 

The velocity field of the SPH particles at t = 7.39 for model B4-1 is shown in Fig. ||, 
where the bar lies along the thick straight line that passes the center. The thin lines represent 
the hypothetical lines connecting the center and the observer. In order to compare with the 
observations (CO and HCN I — v maps), we have constructed theoretical / — v maps by using 
this velocity field of SPH particles. The observer is assumed to move along a circular orbit with 
rotational speed vq = 220 km s _1 at a galactocentric distance of 8.5 kpc. Several / — v maps are 
constructed with different values of which is the angle between the bar and the line-of-sight. A 
representative I — v map at = 50° is shown in Fig. |6|. This figure should be compared with the 
CO and HCN data shown in Fig. H. We pay attention to features such as the forbidden velocity 
components and the velocity gradient produced by the central concentration of SPH particles close 
to the GC. 

The particles with the negative velocity at the positive longitude or the positive velocity at 
the negative longitude are those moving along non-circular orbits. The forbidden components seen 
at I = 1.7° and the Sgr A, C, and E complexes also appear in the SPH simulations. However, 
our SPH results are quite different from the CO map in Fig. [l|b. The CO map simply shows 
much more prominent non-circular components than the simulated results. SPH particles initially 
on circular motions quickly acquire non-circular velocity components in the early phase, but the 
hydrodynamic collisions tend to suppress the random velocity. 
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The steepness of the diagonal feature in Fig. |6| (i.e., the velocity gradient) depends on 
the elongation of the bar. As the bar becomes more elongated, the velocity gradient along the 
longitude of these concentration tends to be steeper. Factors determining the bar strength are the 
mass and the axial ratio of the bar. We assumed the bar mass to be about 5 % of the total mass 
of the Galaxy. Among our models with different elongation, the model B4-1 (b/a = 1/4) produces 
the steepest velocity gradient of central concentration. Since the velocity gradient implied by the 
HCN data is very steep, the elongation of the bar in our galaxy is likely to be rather large. 

The disk component in nearly circular motion appears as a narrow band running from the GC 
to I ~ ±6° with velocity width of |V^5^| ~ — 50 km s _1 . The slope of the disk component in the 
I — v map depends on the extent of how far away we are watching the gaseous material. If clouds 
exist all the way to the solar circle, the slope would be nearly zero. However, if the clouds exist 
only close to the GC, the lower velocity envelope would appear as a straight line with nonzero 
slope. Since the gaseous material exists only within 4 kpc in our simulation, a slope is seen in our 
model I — v map. The actual disk component seen by the CO appears as a strip with nearly zero 
slope. On the other hand, the HCN I — v map (Fig. |l]b) does not show the disk component. This 
means that the HCN clouds are fairly well confined within a small distance from the GC. 

The gas-depletion region within the CR shown in Fig. [| appears between the outer elliptical 
ring and the central concentration in the I — v map of Fig. ||. This structure matches well with 
the one at \l\ ~ 2° — 4° in the CO map (Fig. |l]a). The depletion of particles results from their 
inflow as they lose orbital energy through hydrodynamic collision. We discuss this phenomenon in 
greater detail in §5. 

5. DISCUSSION 

5.1. Collisional and Collisionless Clouds in the GC 

Our models have two important ingredients: a non-axisymmetric potential and hydrodynamic 
effects. The non-circular velocity component is induced by the fluctuating potential of the 
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rotating bar. The hydrodynamic effects generally suppress further developments of non-circular 
motions. In order to distinguish between these two effects, we have made further experiments 
with the "best-fitting" B4-1 model by turning off the pressure force which is responsible for the 
hydrodynamic collisions between the SPH particles. This "collisionless" model resulted in the I — v 
diagram shown in Fig. [5]. Obviously, the rate of inflow of the gas within the CR is much smaller 
than that of the SPH model, and the 'gas-depletion region' is not seen in this figure. The dispersion 
in velocity is much greater than that of the case with collisional effects (i.e., the computation 
including the pressure force). Thus, we clearly see that the large fraction of non-circular orbits 
induced by the potential fluctuation are suppressed by the presence of hydrodynamic collisions. 
Without such collisions, the radial velocity of the gas clouds would have covered a much larger 
area in I — v space. 

We have also performed orbit calculations of single particles to try to understand the 
kinematic structure of the clouds. The hydrodynamic collisions between molecular clouds tend 
to destroy non-closed, self-intersecting orbits. In a rotating, non-axisymmetric potential, there 
are several families of closed, non self-intersecting orbits. Among them, the xl and x2 orbits 
are the most prominent ones (Contopoulos & Mertzanides 1977). These orbits and the resulting 
I — v traces are shown in Fig. ||. The I — v traces are drawn assuming that the angle between 
the bar and the line-of-sight is 45°. The xl orbits shown as solid lines move along the long axis 
and provide a rather large dispersion in the radial velocities when viewed from the position of 
the Sun unless the line-of-sight coincides with the short axis. Binney et al. (1991) have shown 
that the parallelogram structure of CO clouds can be generated by the cusped orbit (the smallest 
non self-intersecting orbit of Contopoulos' sequence xl) in a simple, non-axisymmetric rotating 
potential. However, in our model the xl family of orbits is not present among the SPH particles. 
We have argued that the HCN I — v map is rather similar to our SPH simulations. This can be 
understood if we assume that the SPH particles represent the dense clouds that are formed by 
hydrodynamic collisions between clouds. On the other hand, the collisionless particles can be 
assumed to represent the less dense clouds that can be traced by the CO emission. This implies 
that the xl orbits are easily destroyed by the hydrodynamic collisions. However, the x2 orbits 
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appear to be more prominent in 'central concentration' in the central part. Although the xl 
orbits do not interact among themselves, there must be collisions between particles with xl and 
x2 orbits. 

In order to make a more detailed comparison, we have also made another simulation with 
a mixture of collisionless (about 10,000) and collisional (about 10,000) particles. The resulting 
I — v diagram in Fig. |9| successfully reproduces the observed features such as the parallelogram 
structure, the high-velocity envelopes of |Vl5_r,| ~ 170 — 260 km s _1 , the forbidden velocity 
components of Vlsr ~ 100 km s _1 at the negative longitude, as well as the structures seen in 
HCN map. Some of CO clouds appearing in the I — v map may be the clouds whose orbits are 
perturbed by the rotating bar potential, but not experienced by collisions between them. The fact 
that the HCN I — v map resembles that of the SPH particles while the CO I — v map is more like 
that of the collisionless particles indicates that the dense clouds traced by HCN molecules may be 
formed by hydrodynamic collisions between CO clouds. 

A comparison of the I — v map of the 'mixture' model with observations also indicates that 
the I = 5.5° complex may consist of a spiral arm, the elliptical outer ring, and the collisionless 
components. It is also possible that the I = 3.2° complex may include x2 family components of 
the central concentration as well as all components that the I = 5.5° complex does. Hence we 
suggest that the well-known wide line width of these complexes (e.g., Stark & Bania 1986; Lee 
1996) may be due to a complex superposition of these several kinds of orbit components along the 
line-of-sight. This model also suggests that the clouds in forbidden velocities (|VLsh| ~ — 100 
km s _1 ) in I = 3.2°, I = 1.7°, the Sgr A, C, and E complexes possibly originate from the elliptical 
outer ring components, and/or from infalling spiral components which experienced hydrodynamic 
collisions, and/or from collisionless components with high velocity. 
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5.2. Gas inflow to the Galactic Center 

The inflow toward the GC due to the bar potential is observationally a well-known 
phenomenon. Some of the star burst activity in barred galaxies are interpreted as a consequence 
of such an inflow. The basic mechanism for the inflow is the dissipation of orbital energy of 
clouds through the hydrodynamic collisions between clouds. The collision occurs because of the 
non-circular motions induced by the rotating non-axisymmetric potential of the bar. In our 
Galaxy, the rate of inflow through the ILR has been suggested to be ~ 0.1 M Q yr -1 (Gerhard 
1992). 

For more quantitative analysis, we calculate how the fraction of gaseous mass within r = 0.02, 
0.1, 0.2, and 0.3 changes with time, and the results are shown in Fig. |n]. For the purpose of 
checking numerical error, we have also shown the same quantities for a model without a bar in 
this figure as straight lines. It clearly shows there is no inflow for all the radii throughout the 
evolution in the UB model, whereas a significant amount of gas inflow is seen in the models with 
a bar. Within r = 0.3, the fraction of gaseous mass is nearly constant throughout the evolution, 
but it increases noticeably within r = 0.2. 

It should be noted that the increase of the gas mass within r = 0.1 is much larger than that 
within any other region. Because the inflow within r = 0.2 is small, most of the gas inflow must 
originate from the region between r = 0.1—0.2. On the other hand, the inflow within r = 0.02 
occurs mildly because our resolution is not good enough within this radius. 

The behavior of gas inflow clearly depends on the bar strength in that the amount of gas 
inflow is greater in models with larger bar axial ratio. The inflowing gas mass through r = 0.1 is 
estimated to be about 0.1 — 0.3 M yr _1 for B models. These are very similar to that estimated 
by Gerhard (1992). 
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5.3. Forbidden Components 

We have shown that the forbidden components could be generated by the perturbation of the 
bar potential. It is also possible that the non-circular motions can be generated by other ways such 
as an explosive event from the Galactic Center. However, such non-circular components should 
be easily dissipated by the cloud collisions. Spontaneous generation of non-circular velocities is 
required to explain the forbidden velocity components unless we are observing the event at a 
special epoch. A rotating bar potential appears to be a plausible source for the generation of the 
non-circular velocity field. 

The forbidden velocity component is suppressed by the hydrodynamic collisions. The 
HCN clouds with forbidden velocities are relatively rare. This is consistent with the collisional 
hypothesis for the formation of dense clouds. We also found that the forbidden components 
appear only when lies between about 20° and 50°, even though they should be visible at almost 
all angles of Qf, if they were generated only by the perturbation of the bar potential. Therefore 
comparing the SPH simulations with the HCN data gives better constraints on ©5. 

5.4. Effects of the finite lifetime and the self-gravity of molecular clouds 

We did not take into an account of the finite lifetime of molecular clouds in our calculation. 
The lifetime of the giant molecular clouds is known to be be around (1 — 2) x 10 7 years (Palla 
& Galli 1997), which is somewhat smaller than a orbital time scale (~ 10 8 years) of the bar. 
It is not clear what is the life-time of the CO clouds represented by 'collisionless' particles in 
our simulation. If it is very short, the clouds may not have enough time to acquire non-circular 
velocity component. Then many clouds would not have enough time to experience collisions before 
they turn into stars. However, we still expect that the kinematic structures of the remaining 
clouds can be similar to our simulated results. If the dense molecular clouds are formed by cloud 
collisions, we are selectively seeing the clouds that have undergone collisions. The CO clouds may 
have extensive non-circular motions from the birth because the gaseous material should have been 
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exposed to the potential generated by rotating bar. 

We have also ignored the self-gravity of the clouds. By neglecting the mutual gravitational 
interactions, we have effectively suppressed the possibility of the dynamical instability. The 
lopsided distribution of dense clouds could be due to non-axisymmetric instability, as discussed 
below. We will consider this problem in the future investigation. 

5.5. Unexplained Features 

One of the puzzles regarding the distribution of molecular gas toward the GC region is its 
lopsidedness: more molecular clouds are seen at the positive / than at negative /. The degree of 
lopsidedness is greater when observed with high-density tracers such as HCN and CS than with 
low-density tracer such as CO (Bally et al. 1987; Lee 1996). Clearly, our models are unable to 
reproduce the lopsidedness. An attempt to explain this phenomenon by perspective effect was 
not very successful (Binney 1994). The m = 1 instability may be responsible for the asymmetric 
distribution of the molecular clouds (e.g., Yuan 1998). Further studies including the self-gravity 
of the gaseous component may help to clarify this point. 

The activity-related features such as the radio continuum structures vertical to the Galactic 
plane [see Sofue (1989) for a review] cannot be well explained in the context of our models because 
our simulations are restricted to the disk. The extension to three-dimensional model is necessary. 
Although we have successfully reproduced some important features by using a simple model, many 
of the important features certainly can be explored with more realistic models. 

6. SUMMARY 

We have studied response of molecular clouds to a bar potential by conducting a numerical 
simulation of the Galaxy using a SPH method to try to understand the distribution and kinematic 
structures of the GC molecular clouds. The model Galaxy consists of three axisymmetric 
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components (Logarithmic massive halo, exponential disk, and a compact bulge represented by the 
Plummer model) and one non-axisymmetric bar potential. There are two physical ingredients in 
our models: non-axisymmetric potential and hydrodynamic effects. The bar component tends to 
generate non-circular velocity component while the hydrodynamical effects tend to suppress it. In 
order to understand these effects, we have made several calculations that include the collisional 
SPH particles only, collisionless particles only, and both. Followings are our results that we 
conclude from detailed comparisons of these calculations with HCN (Lee 1996) and CO (Bitran 
1987) observation data. 

(1) The velocity structures of the SPH particles are found to be consistent with ones traced by 
HCN line which are very dense clouds. However, the I — v diagram of CO clouds is quite different 
from that of our SPH results, but more similar to that obtained from the simulations without 
hydrodynamic collisions. The I — v diagram of the mixture of the collisional particles (60%) and 
the collisionless particles (40%) gives better reproduction of kinematic structures of the GC clouds 
traced by CO line, suggesting that the dense clouds seen in the HCN line are formed by collision 
between less dense clouds seen in the CO line. 

(2) Our models give one possible explanation of why the line widths of spectra of the / = 5.5° 
and I = 3.2° complexes in the GC are so broad. The model fitting to the data suggests that the 

/ = 5.5° complex may consist of a spiral arm, an elliptical outer ring, and collisionless components, 
and the I = 3.2° complex may include more x2 family components of central concentration as well 
as all components that the I = 5.5° complex does. Hence the well-known wide line width of these 
complexes may be due to such superposition of many different types of orbit components along 
the line-of-sight. 

(3) Our models also suggest what forbidden components (|Vls.r| ~ — 100 km s _1 ) in 

I = 3.2°, I = 1.7°, the Sgr A, C, and E complexes look like. These can possibly originate from the 
elliptical outer ring components, and/or from the infalling spiral components which experienced 
hydrodynamic collision, and/or from collisionless components with high velocity. We also found 
that forbidden components originating from the hydrodynamic collisions are well reproduced 
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when Ob lies between about 20° ~ 50°. However, the forbidden components originating from only 
the perturbation of the bar potential are seen from all viewing angles. The relative rareness of 
the forbidden velocity components in the HCN clouds, compared with those in CO clouds, may 
indicate that the HCN clouds are suppressed by the hydrodynamic collisions to have lower velocity 
dispersions. 

(4) Some features such as the lopsidedness in distribution of molecular clouds to the positive 
longitude, and activity-related ones as radio continuum structures vertical to the Galactic plane 
are not explained by our simple two-dimensional models. More realistic models which include 
the third dimension, self-gravity of interstellar clouds, and so on will be necessary to study these 
phenomena. 
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Luis Ho for a careful reading of the manuscript. We would like to express a final acknowledgment 
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Table 1. Model parameters 



Model 


a 


b 


Mdisk 


Mb u ig e 


M bar 


R c 


Rd 


r c 


B2-1 


0.22 


0.11 


0.30 


0.076 


0.055 


0.70 


0.32 


0.024 


B2.5-1 


0.22 


0.088 


0.30 


0.073 


0.058 


0.75 


0.34 


0.023 


B3-1 


0.22 


0.073 


0.31 


0.073 


0.053 


0.75 


0.33 


0.023 


B4-1 


0.22 


0.055 


0.33 


0.072 


0.045 


0.75 


0.35 


0.023 


UB 






0.33 


0.079 




0.75 


0.34 


0.023 
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Fig. l.— l-v maps of (a) CO (Bitran 1987) and (b) HCN (Lee 1996) observations. 

Fig. 2. — A rotation curve of the model Galaxy with a bar axial ratio of a : b = 0.22 : 0.055 
(=4:1) and the observed rotation curve (Clemens 1985). 

Fig. 3. — A circular frequency f2(R), a bar pattern speed f2&, and f2(R)±i| of the model Galaxy 
with a bar axial ratio of a : b = 0.22 : 0.055 (= 4 : 1), where k is an epicyclic frequency. 

Fig. 4. — Evolution of particle distribution in the B4-1 model. Typical two spiral arms are visible 
at the early stages of the evolution and remain weak at the later stages. Note that there are 
prominent occurrences of the x2 orbits (central concentration) within r « 0.10, the gas-depletion 
region between r ~ 0.1 — 0.2, and the elliptical outer ring between r « 0.2 — 0.3 throughout the 
entire integration. 

Fig. 5. — Velocity field of the particles at t=7.39 for the B4-1 model. The direction of the major 
axis of the bar is shown as the straight thick line passing through the center, and other lines are 
hypothetical lines toward the assumed directions of the Sun. 

Fig. 6. — Longitude— velocity diagram of SPH particles seen at a viewing angle = 50° at the 
t = 7.39 epoch of the B4-1 model. 

Fig. 7. — The I — v map for collisionless particles. 

Fig. 8. — The xl (solid lines) and x2 (dotted lines) orbits (upper panel) and their I — v traces 
viewed at the solar circle (lower panel). We have assumed that the angle between the line-of sight 
and the bar is 45°. 

Fig. 9. — A run of the B4-1 model with the mixture of SPH particles (60%) and collisionless 
particles (40%) (a viewing angle 6& = 40°). The CO data from Bitran (1987) are superposed with 
contours. 

Fig. 10. — The evolution of gas mass fractions within various radii as function of integration time. 
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